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Abstract 

Using Yang and Yang's particle-hole description, we present a thorough derivation of the ther- 
modynamic Bethe ansatz equations for a general SU(k) fermionic system in one-dimension for 
both the repulsive and attractive regimes under the presence of an external magnetic field. These 
equations are derived from Sutherland's Bethe ansatz equations by using the spin-string hypothe- 
sis. The Bethe ansatz root patterns for the attractive case are discussed in detail. The relationship 
between the various phases of the magnetic phase diagrams and the external magnetic fields is 
given for the attractive case. We also give a quantitative description of the ground state energies 
for both strongly repulsive and strongly attractive regimes. 



PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.IK, 05.30.Fk 



I. INTRODUCTION 



Exactly solvable models of interacting fermions in one-dimension (ID) have attracted 
theoretical interest for more than half a century. Before 1950, it was not clear how to treat 
the Schrodinger equation for a large system of interacting fermions. The first important 
breakthrough was achieved by Tomonaga [1] who showed that fermionic interactions in ID 
can mediate new collective degrees of freedom that are approximately bosonic in nature. 

In 1963, Luttinger Q introduced an exactly solvable many-fermion model in ID which 
consists of two types of particles, one with positive momentum and the other with nega- 
tive momentum. However, Luttinger's model suffers from several flaws which include the 
assumption that the fermions are spinless and massless, and more importantly an improp- 
erly filled negative energy Dirac sea. Mattis and Lieb |3| expanded on Luttinger's work by 
correctly filling the negative energy states with "holes". Before that, Lieb and Liniger {4, 5| 
solved the ID interacting Bose gas with ^-function interactions using Bethe's hypothesis [6|. 
Later McGuire solved the equivalent spin-1/2 fermion problem for the special case where 
all fermions have the same 8 p,„ except one having the oppos.te spin in the repulsive fl and 
attractive [8[ regimes. He showed that in the presence of an attractive potential a bound 
state is formed. Furtherprogress by Lieb and Flicker 9] followed on the two down spin 
problem. In 1967, Yang [10] solved the fermion problem for the most general case where the 
number of spin ups and spin downs are arbitrary by making use of Bethe's hypothesis. At 
the same time, Gaudin [U] solved this problem for the ground state with no polarization. 



Sutherland 



12] then showed that the fermion model with a general SU(k) spin symmetry 



is integrable and the solution is given in terms of k nested Bethe ansatz (BA) equations. 



And in 1970, Takahashi 



13] examined the structure of the bound states in the attractive 



regime with arbitrary spin and derived the ground state energy together with the distribution 
functions of bound states in terms of a set of cou pled integral equations. Using Yang and 
Yang's method for the boson case, Takahashi jl5| and Lai [la, [l^ derived the so-called 
thermodynamic Bethe ansatz (TBA) equations for spin-1/2 fermions in both the repulsive 
and attractive regimes. The spin-string hypothesis describing the excited states of spin 
rapidities was also introduced by both authors. Later on, Schlottmann 18|, [l9| derived 
the TBA equations for SU(k) fermions with repulsive and attractive interactions. See also 
Schlottmann's epic review article on exact results for highly correlated electron systems in 
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ID [20|. 

The TBA equations have been analyzed in several limiting cases, i.e., T — > 0, T — > oo, 
c — > and |c| — )■ oo, where T is the temperature and c is the interaction strength. The 
ground state properties and the elemental charge and spin excitations were also studied for 



some special cases. However, the TBA equations for the attractive regime [19|, |20j are not the 
most convenient for the analysis of phase transitions and thermodynamics. For the attractive 
case, it was shown that the ground state in the absence of symmetry breaking fields consists 
of spin neutral charge bound states of k particles. The repulsive case however consists of 
freely propagating charge states and spin waves with different velocities. The phenomenon of 
spin-charge separation plays a ubiquitous role in the low energy physics of ID systems 21]. 
However, the physics of these models, such as the universal thermodynamics of Tomonaga- 
Luttinger liquids, quantum criticality and the universal nature of contact interaction, are 
largely still hidden in the complexity of the TBA equations. It is thus important to develop 
new methods to extract the physics of ID exactly solved many-body systems in order to 
bring them more closer to experiments. 

Most recently, experimental advances in trapping and cooling atoms to very low tem- 
peratures allow a test of the theoretical predictions made so far. In particular, Liao et al. 



22j experimentally studied spin-1/2 fermions of ultracold 6 Li atoms in a 2D array of ID 



tubes with spin imbalance. The phase diagram was confirmed and it was discovered that a 
large fraction of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)-like phase lies in the trapping 
center accompanied by two wings of a fully paired phase or unpaired phase depending on 
the polarization. This observation verified the theoretical predictions 23^ 26] regarding the 
phase diagram and pairing signature for the ground state of strongly attractive spin-1/2 
fermions in ID. Although the FFLO phase has not yet been observed directly, the exper- 
imental results pave the way to direct observation and characterization of FFLO pairing 
22|. 

In this paper, we derive the TBA equations for a general ID system of fermions with 
SU(k) spin symmetry from Sutherland's BA equations using the same approach as Yang 
and Yang for ID bosons [14J]. Both the repulsive and attractive cases are discussed. We 
also give the exact thermodynamics of the ground state of the attractive and repulsive cases 
in both the strong coupling and weak coupling limits. A general relationship between the 
different magnetic phases and the external magnetic field is discussed for the attractive case. 
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How the external magnetic fields affect the different pairing phases in the attractive regime 
is also addressed. This paper gives a thorough derivation of many results in a recently 



published paper [27| that provides the exact low temperature thermodynamics for strongly 
attractive SU(k) fermions with Zeeman splitting and shows that the system behaves like a 
universal Tomonaga-Luttinger liquid in the gapless phase. 

II. THE MODEL 

The Hamiltonian for the ID iV-body problem is 

H =-^Ld^+9^ £ 6{x i -z j ) + 1 £N i 4 s tfB,B). (1) 

i=l % l<i<j<N i=l 

It describes N fermions of the same mass m confined to a ID system of length L interact- 
ing via a ^-function potential. The first and second terms in the Hamiltonian correspond 
to the kinetic energy and ^-interaction potential respectively. The coupling constant gm 
can be expressed in terms of the interaction strength c = 2/ am as gin = —h 2 c/m where 
am is the effective ID scattering length. For repulsive fermions, c > and for attractive 
fermions, c < 0. We shall assume that the system has periodic boundary conditions i.e., 
ip(xi, . . . , Xi, . . . , xjy) = i/j(xi, . . . , Xi + L, . . . , xjy) where Xi is the position of the i-th particle. 
There will be k possible hyperfine states |1), |2), . . . , \k) that the fermions can occupy. The 
number of fermions in the state \i) is given by iV* while the Zeeman energy e l z is determined 
by the magnetic moment fx l B and the magnetic field B. For brevity, we shall choose the 
dimensionless units H = 2m = 1 for the rest of this paper. 

The wavefunction ip for this Hamiltonian satisfies the symmetry of an irreducible rep- 
resentation = [k nk , (k — l)^" 1 -^ K ) _ _ _ ^2 Ar2_Ar3 ; x^ 1- ^ 2 ], where the iV*s are such that 
iV 1 > A^ 2 > . . . > N K . The Young tableau which corresponds to this irreducible representa- 
tion is given in FIG. [TJ This system has SU(k) spin symmetry and U(l) charge symmetry. 



Sutherland 



12] showed that this problem can be solved by repeatedly using the generalized 



Bethe's hypothesis which was introduced by Yang 10] . To obtain the total momentum 
p = X^fli kj and the energy eigenspectrum E = 5^*=i ^| f° r ^^e system, we need to find a 
set of quasimomenta kj that satisfies the equations 

Mi , _ . (1) . , 

exp(^L) = I P AD 3 = 1>-,N ( 2 ) 

1 = 1 kj A; IC 



U> M) 13) 



N 2 -N 3 



N 1 -N 2 



N*-i-N* 

T 



FIG. 1: The Young tableau which corresponds to the irreducible representation = [k , (k 

^N K ' 1 -N K 2 n2 ~ n3 l^-N 2 -^ 

A^ - + ic' Pr A^ - A« + 2id A$ - A^ - id 



l\ A^ - - S = ~ |i Ai p) - A? - 2 % d M AL P) - A? + x) + id 



1 A f/ " - \^ ! - M A"" - \ 

7= 

a = l,...,M p and p = !,...,« — 1 (3) 



where M := iV, M K = 0, A,- := fc,- and A,- is undefined. The rapidities A^ char- 

J J J J 

acterize the internal spin degrees of freedom. The quantum numbers Mi are defined as 
Mi = Ylj=i Afi +1 and d = c/2. This set of k coupled algebraic equations are called the BA 
equations. 



III. THE ROOT PATTERN 



We shall only consider the strong coupling re gim e where L\c\ ^> 1. For the repulsive case, 



it is easily shown that the roots k~ must be real 



28l | . However the rapidities A^ are allowed 



to take on nonzero imaginary values. It was first suggested by Lai [16| that the rapidities 
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appear as strings in the complex plane of the form 

A(pW = A(p).- + z(n + l-2j)|c / |+0(exp(-5L)) j = l,...,n (4) 

in the thermodynamic limit where N, L — > oo while keeping the ratio N/ L fixed. Here n is 
the length of the string, a labels each individual string and Aq denotes the real part of 
each string. Every string must be symmetric along the real axis, i.e. any complex solution 
is accompanied by its complex conjugate pair. This solution for the rapidities hold as long 
as L\c\ ^> 1. In other words, it holds up to order exp(— SL) where 5 is a positive monotonic 
increasing function of the interaction strength |c|. When the system is in its ground state, 
the rapidities do not form strings. The strings also obey the relation M p = Y^nLi n ^n ) 
where AfI P ' is the number of A^ strings with length n. 

In the attractive regime, it is found that complex string solutions of kj also satisfy the 
BA equations. A system with ^-components of fermions has SU(k) symmetry. The quasi- 
momenta kj may appear as bound states of up to length k. For convenience, we denote the 
number of kj bound states with length 1 < i < k as A^. In the ground state, none of the 
bound states can be broken which means that Ni = N l — N l+1 . A bound state in fc-space 
of length m will take the form 

C j = AL m - 1) +*(m + l-2j)|c , |+0(exp(-<5L)) j = l,...,m (5) 

with real part Ai" 1 ^. In general, a k a bound state of length m will be accompanied by 
a Aa^ bound state of length m — 1, a bound state of length m — 2 and so on until a 
Ai m ^ bound state of length 1. Each accompanying bound state in A^ -space, 

A (2) 

-space, 

-space will share the same real part \ ™ . A graphical depiction is given in 

FIG. 13 

However, strings in A( p )-space do not have to be accompanied by any shorter string. 
Therefore, the difference between our definition of a bound state and a string is that a 
bound state is a string that originates from fc-space, and is accompanied by corresponding 
strings in each subsequent A W -space. On the other hand a string in A W -space characterizes 
the spin excitations that can exist independently in spin sectors. 
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FIG. 2: (Color online) Root pattern for the k = 3 case with N = 23, A 1 = 13, N 2 = 7 and A 3 = 3 
in the ground state (T = 0). Bounds states are formed for charge and spin rapidities and are 
indicated by the dashed boundaries. Accompanying bound states in each subsequent A^ -space 
share the same real part. The maximum length of any bound state in fc-space is k = 3. In an 
excited state (T > 0), strings can also form in A^ -space. 

IV. THE TBA EQUATIONS: REPULSIVE CASE 

The TBA equations which are expressed in the form of dressed energies allow us to pre- 
cisely derive numerous thermodynamic quantities and to analyze the behavior of phase tran- 
sitions at critical points. Explicit expressions for the free energy, grand partition function, 
pressure, chemical potential and so on can be directly obtained from the TBA equations. 
Physically, the TBA equations give us the energy of elementary excitations above the ground 
state. There are several steps to take in order to derive these equations. We will give an 
outline of each step in deriving these equations for the repulsive case, all of which follow 



from Yang and Yang's pioneering work 



The string solution for A^ gives us a different form of the BA equations when substituted 
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into the original equations ([2]) and <^fy. After lengthy calculations, we obtain 



oo M™ , .(l), n 



^) - n n (6) 

n =l c=i -^a 

iV .(l),n , ■ ,00 Mi 2 ' oo M™ 

n i." " c , n n ( a ?'" - = - n n £ »»> w 1 '- - m 

i =1 A a K\ inc m=1 p =1 m=1 ( g =1 

oo M< p) 



n n ( A ^ ),n - A ? ),m ) 

m=l /3=1 

oo A4T 1 ' oo Mi p+1) 

n n f «» ( a « ),n - A r i),m ) n n f - ( a « ),n - A ? +i),m ) ( § ) 



m=l /3=1 m=l /9=1 

oo Mi K ~ 2) oo M* K_1) 



n n f - ( A « K_i) * n - A r 2) ' m ) = - n n ^ ( A r i),n - A r i),m ) . (9) 

m=l /3=1 m=l /3=1 

The functions E nm (x) and F nm (x) are 

, i e| n _ m |(a;)e^_ m | +2 (x) . . . e,^ +m _ 2 (x)e )1+m (x), for n 7^ m; 

&nm{X) — \ (,-1-UJ 

el(x)el(x) . . . el n _ 2 (x)e2n(x), for n = m. 

ei n _ m | + i(x)e| n _ m | +3 (x) . . .e n+m _ 3 (a;)e n+m _i(x), for n 7^ m; 

f»W = S ( U ) 

ei(x)e 3 (x) . . . e 2 „_ 3 (a;)e2n-i(a;), for n = m. 



where 



Taking the logarithm of each equation yields 



*) = — 7- (I 2 ) 



E * f^f* ) +E E r- PHrH - £ e„ m ( 

Z=l \ / m=l /3=1 \ / m=l /3=1 \ 

(14) 



00 M% 1} /^(p),n _ ^(p-l),m\ 00 a4? +1) _ ^(p+l),m 



^ ^ r nm I ^ ) + rm7 ' 1 r ' 

m=l /3=1 V / 771=1 (9=1 







2*JW.» + £ £ 6 nm ( - - g ) (15) 

771=1 /3 = 1 
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x> Mfn 2 ' /A( K_ l). n A( K_2 )i m \ oo Afm ^ / A («— l),n a(k-1),?TI n 



EE r - - — ^ — = Mr 1 * +E E e » 

m=l /3=1 \ / m=l /3=1 



(16) 

where Jj and Ja n are odd or half-odd integers depending on the quantum numbers. The 
functions Q„ m (x) and T„ m (x) are 



q / \ I y |n— m| y VJ' 1- m\+z j \n+m—z/ \n+mj • ' > (17) 

20(§)+20(f) + ... + 20(^)+0(f;), forn = m. 

+0(f) + ... + e(-2_)+6l(-^) ) forn = m. 

where 

6{x) = 2tan _1 (x). (19) 

Let us now introduce continuous monotonic increasing functions f(k) = 27rIi/L and 
9n\k) = 2irJa^' n I L. Denote p(k) and p h (k) as the densities of "particles" and "holes" 
in /c-space. Similarly, denote and <Jn^' h as the densities of "particles" and "holes" for 
strings with length n in the p-th rapidity space. We shall express equations (fT3j) to fTT6]) in 
the continuum limit as 

/(*) = k + jr[6 f^f] a£\k')dk> (20) 

n=l J \ nc / 



• \ / m=l ^ ' 



E l®™{ k —f)^{k')dk> (21) 



m=l 



<?i p) W =jt[r nm (LJ£\ at 1] (k')dk' + f; /" r nm f*-^) a% +1 \k')dk' 

"E hn m ( k -^-]at\k')dk' (22) 



m=l 



^H*) = E / r - ( - E /q- (^) 'St^wm- (23) 

m=l J \ S m =\ J \ / 
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Using the relations ±f(k) = 2n(p(k) + p h (k)) and ±$\k) = 2n(a^\k) + a { ^ h (k)) in 
the thermodynamic limit, we then have 



oo 



P (k)+p h (k) = ±- + J2*n*^\k) (24) 

n=l 

oo oo 

° { n\k) + ^\k) = a n * p(k) + J2S nm * <r£\k) - Turn * *£> (k) (25) 

m=l m=l 

oo oo oo 

+ = S nm * ^(k) + Sum * CJ% +1 \k) - £ T nm * ^(fc) (26) 

m=l m=l m=l 

oo oo 

m=l m=l 

where the convolution integral / * g(t) = J f(r)g(t — t)cLt is used. The functions T nm (x) 
and S nm (x) are 

I a| m _ n |(x) + 2a| m _„| +2 (x) + . . . + 2a m+n _ 2 (x) + a m+n (x), for n ^ m; 
T nm {x) = < (28) 
I 2a 2 {x) + 2a 4 (x) + . . . + 2a 2n _ 2 (x) + a 2n (x), for n = m. 

J a| m _„| + i(x) + a| m _„| +3 (a;) + . . . + a m+n _ 3 (a;) + a m+ „_i(a;), for n 7^ m; 
I ai(x) + a 3 (x) + . . . + a 2n _ 3 (x) + a 2n -i(x) : for n = m. 

where 

, . 1 nlc'l 
a " (x) = + (30) 

The Gibbs free energy per unit length is given by Q/ L = Ej L — pN/L — TS/L + E z /L 
where the first term corresponds to the ground state energy, the second term corresponds 
to the addition and extraction of particles from the system, the third term represents the 
entropy of the system and the last term is associated with the Zeeman energy per unit length 
of the system, p is the chemical potential, T is the temperature and S/L is the entropy per 
unit length. For each infinitesimal interval dk, the energy is degenerate. The total number 
of microstates with this energy degeneracy is given by 

(Lp(k)dk + Lp h (k)dk)\ ^ ^ (Latt\k)dk + La^' h (k)dk)\ 
(Lp(k)dk)\(Lp*(k)dk)\ * 1111 (La£\k)dh)\(LvP*(k)dk)\' [ } 
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Using Sterling's approximation mm! = mlnm-m when m >> 1, the entropy per unit length 
of the system follows as 

| = / p")ln(P+ (>")- f>hu>- p"hip"\(W 



= J {(p + p h )\n(p + p h )-p\np-p h \np h } 

OO K — 1 n 

+ E E / ( ( ff n ) + + ''"J - ^ 111 ^ - ^ ^ 111 ^ } dk ( 32 ) 

n=l p =l ^ 



where the entropy is defined as 5 = f IncW. 
The Zeeman energy can be expressed as 

K-l 

E z = ~Ys H p NP 



P =i 

K-l 

-ifxJV + Y,(H P ~ H P+ i)M p . (33) 
p=i 



(34) 



The following identities are required for further derivation 

j = [ p(k)dk, | = I k 2 p(k)dk, ^ = E / na^(k)dk. 

J J n=l 

In the thermodynamic limit, the Zeeman energy can be written as 

771 p ^ 1 p 

-f = -Hxl p{k)dk + E^ " ^ih-i) / na%\k)dk. (35) 

^ p=l n =l J 

Minimizing the Gibbs free energy per unit length when the system is in equilibrium i.e., 
dVt = finally gives the TBA equations 

OO 

e(k) = k 2 - fj, - H x - T a n * In (l + e ~«" )(fc)/T ) (36) 

71=1 

OO 

e\k) = n(H x - H 2 ) - Ta n * In (l + e-^/ T ) + T £ T nm * In (l + e~^^ T ) 

m=l 

oo 

-T X> m * In (l + e-^)/^ (37) 



m=l 



OO OO 

#>(*) = n(ff p - //p +1 ) + T £ T„ m * In (l + e"^^) - T £ S nm * In (l + e^"^) 

m=l m=l 

oo 

-Tj2S nm *ln(l + e-^ +1) ^ T ) (38) 



m=l 
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oo oo 

e~ 1} (k) = nH^t + T £ T nm * In (l + e^^) - T £ S nm * In (l + e^"^^) 

m=l m=l 

(39) 

where we have defined p h (k)/p(k) = exp(e(k)/T) and (Jn^ h {k)/(Jn^ = exp(C,n\k)/T). 
The pressure per unit length of the system is 

an T 



P dL 2tt 



J In (1 + e" £(fe)/T ) dife. (40) 



V. THE TBA EQUATIONS: ATTRACTIVE CASE 



After substituting the complex solutions for the bound states and strings into the original 
BA equations (|2j) and (|3]), we obtain 



exp(zmAf)L) = (1)^ 1 II 1111 



' Vg Xf ] -A{ 9) -i(q + m-2p)\J\ 
P =t q=P t=i X i r ) -^ ) + ^ + m-2 P )\c'\ 

^ ' J^il\ ^ - X\ 9) + i(q - m)\d\ 
oo M n A (m) _ A (m),n _ m|c/| 



X 



TTTT^h n — for m=i,...,K (41) 



OO M„ 



(p),m\ 



m=l 3=1 



»(p),n \(p) I • I /i 00 M ™ 

n X'V n n f ->w - A r i,,m ) 

l=l i\ a a 1 m\c | m=1 p =1 

oo M m 

x nn F -( A « ),n - A ? +i),m ) f ° r p=i,...,*-i. (42) 

m=l 0=1 

Taking the logarithm of each equation yields 

m-X k N q / ,( m ) _ ,(q) \ K AT, /\( m )_\(?)' 

J ^^tT \(q + m-2p)\c'\) q ^ + ^ \{q-m)\d\ 
oo M n / \(m) A (m),n\ 

+ EE/^M (43) 

n=l a=l \ 11 / 
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Nr > / ' (p),n \ (p) ' 



E* 



n c 



oo M m /«(p),n A (p-l),n\ oo M m / a (p),n a (p+l),n ' 



- v v r „ m -p ! vv r „„, , 

m=l /3=1 \ 11 / m=l (3=1 V ' 

oo M m / *(p),n _ A(p),m\ 

+ EE e - - u. g • ( 44 ) 

m=l /3=1 V 11 / 

We introduce the continuous monotonic increasing functions h m (k) = 2txK^ jL and 
jn\k) = 2nLa '' n /L. Denote p m {k) and p^k) as the densities of "particles" and "holes" 
for the bound states with length m. Similarly, denote and <Jn^' h as the densities of 
"particles" and "holes" for strings with length n in the p-th rapidity space. In the continuum 
limit, equations (143]) and (jH]) become 

h m (k) = k m-2p)\A ) " ,ik ' )dk ' 

- ± J" ((j^f) *<*0* - E/« (^) "SW («) 

g =7n+l* / \« /II/ n= i •/ \ 11/ 

= f 6 Pp (k')dk' + f; |r nm f*^) a^~ 1 \k!)dk' 

\ 11/ m=l \ I I / 

00 " / T T / \ OO / U 1J\ 



+E/ ^(^e^w^-E/ e -(Vr) 



(k )dk . (46) 



From the relations ^M*0 = 27r(p m (&) + ^ (A;)) and ^jn\k) = 2n(a { ^(k) + ai p),/l (A;)) 
in the thermodynamic limit, we obtain 

m— 1 K K oo 

Pm(fc) + Pm(&) = ^ - ^^a ?+m _ 2p *p ? (A;) - a q _ m * p q (k)-J2 a n*vt\ k ) ( 47 ) 

p=l <?=p q=m+l n=l 

oo oo oo 

+ 4 P) '"(^) = «n * Pp(fc) + S nm * ^ ' W + E ^ * ^ W " E T ™ * ^ (*)" 

m=l m=l m=l 

(48) 

The Gibbs free energy for the attractive case has the same expression as the repulsive 
case. However, the expressions for each term in the Gibbs free energy is different from the 
repulsive case. The total number of microstates with the same energy degeneracy in the 
attractive case is given by 

dw _ TT ( L Pm(k)dk + Lp h m (k)dk)\ K fl (Lalr\k)dk + La^ h (k)dk)\ 

' 11 (L Pm (k)dk)\(Lft(k)dk)\ ^l±{La^\k)dk)\{La^(k)dk^ 1 } 
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The entropy per unit length of the system is 
S ~ f 

T = ^2 {(Pm + fL) ln (Pm + PL) - Pm hi Prn ~ pt hi P„} dk 

m=l 

oo re— 1 „ 

+ EE / {(^+^M^+aW fc )-(7WlnaW-aW*lna« fc }dfc. (50) 

n=l i=l ^ 

The ground state energy per unit length was given by Takahashi [3] as 

f = E/ ( rf " ^^Icf ) M*)*. (51) 

m=l ' 

It can be easily derived by taking the discrete sum i£ = ^ ■ fc| with the roots given in 
equation (J5J) and then extending it to the continuum limit. Therefore the Zeeman energy 
per unit length is 

re-l 

L 



Ez rr 
/ j Timlin 



m=l 

re— 1 « re— 1 oo « 

= -J2 H mj Pm(k)dk + J2( 2H ™ ~ H ™~ i " H ™+^ / nat\k)dk (52) 

m=l ' m=l n=l 

where if K = because we only need k — 1 independent parameters to describe the relative 
distances between the energy levels of different fermionic species due to Zeeman splitting. 
Here we shall denote n m = N m /L for brevity. 

Minimizing the Gibbs free energy with respect to deviations in the various densities yields 
a set of coupled integral equations. On introducing the dressed energy terms exp(e m (/c) /T) = 
Pm(k)/Pm(k) and exp(£n\k)/T) = <Jn^' h (k) / <7n\k) , we arrive at the TBA equations for 
attractive fermions with arbitrary spin 

/ 2 i \ m—1 re 

e m (k) = mk 2 - mfi -H m - ml ™ ~ l) \cf + T V V a q+m . 2p * In (l + e~^ T ) 

K OO 

+T a «-™ * ^ (1 + e~^ T ) - T £ a n * In (l + e-^«/ T ) (53) 

(j=m+l n=l 



oo 

= n(2H p - Hp-i - H p+l ) + Ta n * In (l + e~ e ^ T ) + T ]T T nm * In (l + e -€^(*)/ r ) 

m=l 

oo oo 

-T £ S nm * In (l + e^WA 1 ) - T £ S nm * In (l + e^ +1) «/ r ) . (54) 



m=l m=l 
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Take note that from the definition o4 = <Jn = given earlier, (k) is undefined. Here 
£m{k) with 1 < m < k are the dressed energies for the bound states of length m e.g. £\{k) 
is for unpaired fermions, £ 2 {k) is for pairs, £3 (k) is for trions and so on. The Fermi level 
is at e m (Q m ) = which implies that the bound states of length m are only occupied with 
fermions having quasimomenta — Q m < k < Q m . There is an equivalent description for the 
band fillings of the strings. 

The pressure per unit length for the system is 



VI. THE GROUND STATE: STRONG COUPLING LIMIT 

A. Repulsive case 

In this section we present the thermodynamic properties of the system in the ground state 
where T — > and also in the antiferromagnetic case where H p = for every p = 1,2,..., /c — 1. 
In this regime, there are no string solutions for each level of rapidity because strings only 
exist in excited states. Hence the TBA equations simplify to 



In the regime c > 1, we can estimate a x *e(k) « —2irPai(k) from equation fj40|) . Taking 
the Fourier transform of equations (IBTj) to ( J59l) yields the difference equations 




K 



(55) 



e(k) 

Z (p) (k) 
&- 1] (k) 



k 2 - fj, + <Zi * £ (1) (*0 

a x * e(Jfe) - a 2 * ^\k) + a x * ^ 2 \k) 

-a 2 * £ ip) (k) + a x * ^-^(k) + fll * £ {p+1) (k) 

-a 2 *e (re " 1) w+«i*e ( ' t " 2) w- 



(56) 
(57) 
(58) 
(59) 



= S(oj) [-2t;P + $ 2 \u) 



(60) 



(61) 
(62) 



where 



S(u) 



1 



(63) 



2 cosh \ojc'\ 
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The general solution to this set of difference equations is 

= 2 ^ Psinh ((^-:P)M) 

siiih(/s:|a;c / |) 

The inverse Fourier transform for this function is 



7T 



P 



(64) 



(65) 



Kd \ cos(^) - cosh(4 
The convolution integral in equation (1561) can be decoupled when c ^> 1 where it becomes 



Using Parseval's theorem, 



we obtain 



nP f 1 nd sin(^) 



Kd J 7T (nc') 2 + A; 2 cos(f ) - cosh(g 



-<i/c. 



/(t)p(t)tft = — / f(co)g(-co)du 



(66) 



(67) 



a x *^ (1) (A;) 



-P 



-2P 



e -\uc'\ sinh((«- 1)1^1) ^ 



sinh(K|a;c / | 

oo e ( K -2)\wc'\ _ e -K\uc'\ 
pK,\uic'\ p — n\u)c'\ 



du 



—,[C ( - 



(68) 



where C = 0.577... is the Euler-Mascheroni constant and *ff(x) is the Psi (Digamma) 
function. The values of for k = 1,2,3 are = -C, #(§) = -C - 2 In 2 and 

\&(|) = — C — | In 3 — Therefore equation ( |56l) becomes 



(69) 



Using the conditions e(±Q) = 0, —2ttP = J_ Q e(k)dk and n = dP/dfi followed by 
iteration to keep terms up to order 1/c gives the thermodynamic expressions 



Q irn 



1 + — C + tf - 



2 2 

/i ?s 7r n 



1 + 



P 
P 



2 2 2 

3 

-vrV 
3 



16 

3ft7 
6 



c + m ( - 



1 + — ( C + tf ( - 
1 + —[C + V[ - 

K7 



(70) 
(71) 
(72) 
(73) 



where 7 = c/n. 
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B. Attractive case 



Having derived the expression for the densities of bound states in equation (|47j) . we can 
derive the ground state energy from equation (1511) . As with the repulsive case, there are no 
string solutions in the ground state. In the strong coupling regime |c| 3> 1, equation (HTj) 
simplifies to 



m n q (q + m — 2p)\c' 

2tt ~ ^ ^ ~k (q + m - 2p) 2 \d\ 2 - 

p=l g=p v 

An expression for the Fermi points Q m can be derived by first evaluating the relation 



Pm{ > 2ir 2-^2-^ vr (q + m-2p) 2 \d\ 2 + k 2 ^ w („ _ m \%\,i\% + fc2 ' 1 ' ' ; 

p=l <j=p q=m+l 



n g (q — m) 






7r (g — m) 2 


c'| 


2 + P 



n 



7)1 



|_r p m (k)dk which gives 

TnQm 2r?,q ^_ 1 / Q m \ y-\ 2n q ^_ 1 / Q m 

" ^ V(g + m-2p)|cV V ^ V(g-m)|c'| 

_ mQ m / y^^A 2n q y^ 2n q \ 

and then rearranging to obtain 

= iirirn L + y^ y^ 2n g + -A 2n g \ + Q / J_\ 

m \ ^-^ m(q + m — 2p)\d\ m(q — m)\d\ ) Vlc'l 2 / 

\ p=l q=p q=m+l v / 

The ground state energy per unit length fl5Tl) is then given by 

73 [ Qm dk V f Qm mn q (q + m - 2p)\c'\k 2 

T ~ J 2tt ~ ^ 2-*> 2-" J o vr (o + m - 2p) 2 \d\ 2 + k 2 

-E E ^^^-t/^iw 

2mn q Q z m -A A 2mn q Q z m 

^ 3vr ^-J ^-f ^ 3vr(g + m- 2p)|c'| ^ ^ , 37r(g-m)|c'| 

m=l m=l p=l 9=P m=l g=m+l v /i i 



Em{m 2 — 1) /|2 



3 

772=1 



' 3m ^— ' ^— ' ' 3m 2 (q + m — 2p)\c\ ' ' 3m 2 (q — m)\c\ 

m=l m=l p=l <?=p m=l g=m+l v 1 

_^™<™i z i) nm|c|2 + m (79) 

m=l \ I I / 

This expression does not include the Zeeman energy, which is equal to — Ylm=i n mH m . 
The actual ground state energy in the presence of an external magnetic field must include 

17 



this term. One can easily derive the ground state energy up to arbitrary orders in 1/c by 
including higher order contributions to the Taylor expansions of the functions we applied 
them to. For better accuracy, we derived the ground state energy that includes terms up to 
order 1/c 2 in the more compact form 

E vr 2 < / 2A m 3A 2 \ ^ 

m=l v 11/ rra=l 

where 

m—l k . k . 

Am = yy _^Jih (si) 

z — ' m[q + m — 2p) z — ' m(q — m) 

p=l q=p ^ ?! q=m +l v ^ 1 

and e m is the binding energy for a bound state with length m i.e., 

m(m 2 — 1 



c\\ (82) 
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Looking back at the first TBA equation fl53|) for the attractive case, we can denote 
fi m = /i + Sm. -f Sm. as the effective chemical potentials for the bound states of length m. An 
expression for the effective chemical potentials can be derived as 



I'c 




1 + ^ + Ml") + V" 27r 2 n^ <9A m + 27T 2 n^ m M r 



a 2 \ |c| |c| 2 7 ^— ' 3ma|c| dn a ^— ' malcl 2 <9n„ 

\ ii ii/ rn=l m=l 1 1 

= ^l( 1 + ^ + ^) + ^ + !^. (83) 

a V l c l \ c \ J \ c \ \ c \ 

Here we used the vector notation A = (A\, A 2 , . . . , A K ), B a = (B^, B 2 , . . . , B*) and / = 
(1,1,..., 1) where 

2n 2 nt dA m 



Bl 



3ma dn a 

l-K 2 n z m ( Q(a — m — I) ^-4 Q(a — j) 



+ y ^--» (84) 



3m 2 a y (a — m) ^-^ (a + m — 2j) 

for m, a = 1,2,..., re. The function 0(x) is the Heaviside step function with properties 
G(x) = when x < and 0(x) = 1 when x > 0. For the special cases of SU(3) and <S77(4) 



fermions, see refs. 



29 



30]. 



Zeeman splitting can be characterized by the parameters e™ or H m . e™ is the Zeeman 
energy level for the species of fermions in state \m). H m on the other hand parameterizes 
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the Zeeman energy level for bound states with length m. Both sets of parameters are related 
by the expression 

K K — 1 

~~ ~~ m+1 ). (85) 



J2^z nm = -J2 H ^ ur - " 

m=l m=l 

A consistent solution to this relation for all n m is e\ 



z — —Hi, e z = Hk-i and = H m _i — H m 
for 2 < m < re — 1. If we denote the difference between the energy levels of fermions in state 



\m + 1) and \m) as Am+y 
I A 2il \ 



m+l 



~-Z i 



we obtain the matrix relation 



A 3)2 
A 4l3 



V A Kft — i y 



/ 2 -1 
-1 2 -1 
-1 2 



W #1 \ 

#3 



-1 



V 



(86) 



2 -1 
-1 2 / 

where the blank entries in the upper and lower triangular sections of the (re — 1) x (re — 1) 
matrix are defined as zero. 

I — I 

A useful relation between H m and the effective chemical potentials is [271 ] 



m{fi m - fi K ) H e m 



(87) 



from the original definition of the effective chemical potential. The external fields H m can 
be tuned experimentally to drive the system between different phases where bound states of 
different lengths exist. In the special case called pure Zeeman splitting where A m+ i /m = A 
for all m, the system has three distinct magnetic phases. The first phase consists of only re- 
bound fermions when Hi < H° . The second phase contains a mixture of re-bound fermions 
and unbound fermions when H^ 1 < Hi < H± 2 . And the third phase is made up of only 
unbound fermions when Hi > H^ 2 . 

The critical external fields H^ 1 and H^ 2 can be evaluated from equation (IHTl) by taking 
m — 1 and using the expressions for the effective chemical potentials from equation ()83l) . 
Doing so gives 



Hi 



n 



re 



12 



| 7 | 2 + 7r 2 (m z Y 1 + 



32m 2 



7T 



- — (1 -m z ? I 1 



re(re — 1)|7| 3re(re — 1)|7| 

32m 2 16(1 — m z 



3re(re-l)| 7 | 3re(re-l)|7| 



3re 2 



K-l 1 ' 



(88) 
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where m z is the spin normalized magnetization per particle density i.e., m z = n ^ 2 _^ M t 2 rue . 
This means that while M* Tne depicts the true magnetization of the system, m z normalizes it 
and only takes on values from to 1. The factor corresponds to the species of fermion 
that has the highest hyperfine spin because the unbound phase is made up of these fermions 
only. 

The critical field iff 1 corresponds to m z = 0, while iff 2 corresponds to m z = 1. Substi- 
tuting these values for m z into the equation for Hi yields the general results 



cl 



iff 



c2 



n 



n 



K 



12 



71 



+ 



16 



K 1 



| 7 | 2 + 7T 2 1 



3k(k — 1) 7 3k 2 | 7 



re— 1 



(89) 



(90) 



12 1 " V 3k(«-1)|7| 

The system has a linear field-dependent magnetization near the critical points. For a 
field slightly above iff 1 , the magnetization is given by 



<Hi-H?) 1 + 



re— 1 



27T 2 n 2 v x L \ k(k — 1)| 7 | K 2 \l\~[j 
On the other hand, for a field that is slightly below iff 2 , the magnetization is given by 



(91) 



m 



2n 2 n 2 



k(k- 1)| 7 | 



(92) 



VII. CONCLUSION 

We have presented a thorough derivation of the TBA equations for a system of ID multi- 
component 5-function interacting fermions in the presence of external magnetic fields. The 
key results, in terms of which the thermodynamic properties are obtained, are equations 
for the repulsive case and equations (j55|) - (j53|) for the attractive case. The form of 

20], 



our TBA equations differs from those derived by Schlottmann 



3-20 



possibly interchangeable. To see how this can be done for the SU{2) case, the reader is 



referred to Takahashi's book [28]. The nature of charge bound states describing different 
sizes of atomic molecules was studied in terms of BA root patterns in the attractive regime. 
Quantum phase diagrams and quantum phase transitions were analytically studied from 
the dressed energy formalism. We also presented the ground state energies for the strongly 
repulsive and strongly attractive regimes. We found that all phase transitions for ID 8- 
function attractive fermions are of second order with a linear field-dependent magnetization 



20 



in the vicinities of critical points, rather than a square-root field- dependent magnetization 



18 



-|20|. The 



model 



31 



32|. 



inear field-dependence is as found for the SU(2) case 



231 ] and for the Hubbard 



For the SU(2) case, the TBA equations provide a comprehensive understanding of FFLO 



pairing and finite temperature thermodynamics of Tomonaga-Luttinger liquids [23j, |33J . The 
key features of phase diagrams and low temperature density profiles of trapped ID spin-1/2 
fermions were experimentally confirmed by matching theoretical predictions from the TBA 
equations, see Liao et al. 22| and references therein. For the SU(k) case, these TBA 
equations were used to study the universal thermodynamics through the derivation of the 
equations of state |27j. The results presented in this paper provide the setting for further 
study of quantum critical behavior in ID interacting Fermi gases, where the exact BA 
solutions provide insight into the physical origins of quantum criticality. 
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